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Abstract 



We present a microscopic model for a'-NaV205. Using an extended Hub- 
bard model for the vanadium layers we derive an effective low-energy model 
consisting of pseudospin Ising chains and Heisenberg chains coupled to each 
other. We find a "spin-Peierls-Ising" phase transition which causes charge 
ordering on every second ladder and superexchange alternation on the other 
ladders. This transition can be identified with the first transition of the two 
closeby transitions observed in experiment. Due to charge ordering the effec- 
tive coupling between the lattice and the superexchange is enhanced. This is 
demonstrated within a Slater-Koster approximation. It leads to a second in- 
stability with superexchange alternation on the charge-ordered ladders due to 
an alternating shift of the O sites on the rungs of that ladder. We can explain 
within our model the observed spin gap, the anomalous BCS ratio, and the 
anomalous shift of the critical temperature of the first transition in a mag- 
netic field. To test the calculated superstructure we determine the low-energy 
magnon dispersion and find agreement with experiment. 
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I. INTRODUCTION 



The layered oxide a'-NaV205 has attracted great interest since 1996, when Isobe and 
Ueda reported a phase transition at T = 34 K with a spin-Peierls hke spin gap formation [|l| . 
At low temperatures the spin gap has a size of about A ^ 100 K which yields a BCS 

ratio 2A//cbT(7 ~ 6, much higher than for other organic or inorganic spin-Peierls materials, 
for which it lies around the canonical BCS-value of 3.5. 

Furthermore, experiments have shown that there are actually two transitions, which lie 
very close to each other 0,^. Both are of second order. The first one at Tci ~ 34 K is 
accompanied by a logarithmic peak in the specific heat ||^,[^ while the second at Tci — Tc2 ~ 
0.3 K is of mean-field character evident from a jump in the specific heat. NMR measurements 
suggest that the first transition leads to charge ordering while the second one opens a spin 
gap §. 

Measurements of the critical exponents yield jSs ~ 0.15. . .0.2 for the critical exponent 
of the lattice distortion 6 and /5a ~ 0.34 for the critical exponent of the spin gap A 



From these values the existence of two transition can also be inferred indirectly. Close to the 
critical point the spin gap A due to a lattice distortion S is expected to obey A oc 5^^'^ []12 



corresponding to /3a = 3/4/^5. This relation is not fulfilled in a'-NaV205, indicating the 
existence of two separate transitions. 

For T > Tci a'-NaV205 has equivalent V sites |1T3H15[] implying valence 4.5+. For 
T < Tc2 ^^V-NMR measurements so far show only two inequivalent sites []TB| while x-ray 
structure determination reports three inequivalent sites ||T7|-[19|. 



An interesting experimental observation is the magnetic field dependence of Tci- It is 
only about 25% of the value expected for a spin-Peierls transition [^,20,21| but on the other 
hand, it is much higher than what would be generally expected for a structural transition p2[ . 



In this article we explain a number of these above mentioned features by an analysis 
of a microscopic model. We find that the phase transition at Tci can be regarded as a 
combination of charge-ordering on every second ladder and superexchange alternation on 
the other ladders, a "spin-Peierls-Ising" transition. This also explains the anomalous shift 
of Tci in a magnetic field. The phase transition at Tc2 can be regarded as a spin-Peierls 
transition on the charge-ordered ladders. It is driven by the charge ordering. The lattice 
distortion accompanying the charge ordering increases the coupling constant of a lattice 
mode invoking superexchange alternation on the charge-ordered ladders. With increasing 
charge ordering for decreasing temperatures the coupling constant increases until the second 
phase transition takes place at Tc2- This second phase transition opens a spin gap in 
agreement with experimental observations. Charge ordering is not yet complete at Tc2- 
The coupling constant therefore continues to increase and due to this the low-temperature 
spin gap is larger than what would be expected from Tc2, i-e., the BCS ratio is enhanced. 

We test our low-temperature structure against x-ray structure determination and find it 
to agree within experimental resolution. Furthermore we calculate the low-energy magnon 
dispersion and find it to agree nicely with experimental results from inelastic neutron scat- 
tering 

This article is organised as follows. In the next section we consider a microscopic model 
which we project onto a spin-pseudospin model describing low-energy charge and spin degrees 
of freedom similar to [^^-p5|. We find that a "zig-zag" charge ordering as observed in 



2 



experiment requires the inclusion of spin and lattice degrees of freedom, not only the charge 
degrees of freedom. We therefore extend the spin-pseudospin Hamiltonian to include lattice 
degrees of freedom. 

In section |II| we analyse this Hamiltonian using RPA. We find two phase transitions oc- 
curing one after the other, one being a combination of charge ordering and spin-dimerisation 
on different ladders, the other being a pure spin-Peierls transition. 

In section ^ we compare the calculated low temperature structure with experiment. 
As regards x-ray structure determination we argue that within experimental resolution the- 
ory and experiment agree. With the calculated low-temperature structure we explain the 
magnon dispersion observed in experiment. Finally we summarize and discuss our results. 



II. DERIVATION OF THE MODEL HAMILTONIAN 



To construct a model Hamiltonian we note that according to LDA+U calculations 
the orbitals around the Fermi level are mostly of vanadium dxy character. They are separated 
both from the lower lying oxygen p orbitals as well as from the remaining vanadium rf-orbitals 
and the sodium 3s orbital which lie energetically higher. This means that we have to consider 
only the quarter filled d^y orbitals. For a discussion of charge ordering we restrict ourselves 
to a 2D model neglecting hopping matrix elements and Coulomb interactions between the 



vanadium layers. The former have been shown to be negligible both by LDA-I-U p6[ and 



Slater-Koster type calculations |T9[. The Coulomb interactions are expected to be small due 
to the large distance between the V sites of neighbouring layers. They are also screened 
by intermediate O ions. Therefore our inital Hamiltonian represents an extended one-band 
Hubbard model at quarter filling, taking on-site and intersite Coulomb interactions into 
account: 

H = tR (al^ttj^ + h.c. j + ^ (alaO-jo- + h-c j 

{id)iL <.i,j:>D 
+ VRUirij + ^ VLUirij 

+ ^ ViLnifij + ^Uni^^ni^. (1) 

Here (, )/l denotes pairs of nearest neighbour vanadium c?a;y-orbitals while <^,^r and 
■C, denote pairs of next nearest neighbour vanadium (ia;j^-orbitals along the rung (R) or 
the leg (L). The first four terms describe the effective hopping of the d-electrons between V 
sites i and j. This hopping can take place both via oxygen orbitals and sodium orbitals. The 
last four terms describe the intersite and on-site Coulomb repulsion. These terms connect 
sites as shown in Fig. 0. 

At T > Tci the system is at quarter filling and the V sites are all equivalent [jl3 



Therefore there is an average of one electron per rung. The on-site |ll5|j27|| and the intersite 



Coulomb repulsion create a charge transfer gap, causing Q;'-NaV205 to be an insulator. 



3 



Therefore hopping between rungs takes place only virtually. This enables us to use an 
effective Hamilton operator Hpp which acts on the Hilbert subspace V of all states with one 
electron on each rung. Let Q be the Hilbert subspace complementary to V and P be the 
projector onto V and Q the projector onto Q. An eigenstate iV'pp) of Hpp satisfies: 

Hpp\ijfp) = {PHP - PHQ {QHQ - ^o)"' QHP) li^^p) 

= E,\4l). (2) 

To calculate Hpp we need to know the energies Eq which can be either determined self- 
consistently or, as for a Schrieffer- Wolff transformation, set equal to the eigenenergy Eq^ 
of PHP. In the following we are interested in Hpp only to order O {t'^/E), i.e., to second 
order in the hopping. In this case we can neglect the hopping terms in [QHQ — Eq), since 
both PHQ and QHP scale with t. 

Within the subspace V we use the operators a\^„, aiaa for the electrons in the atomic 
orbitals. Here i denotes the rung, the pseudospin variable a = ±i describes whether the 
electron occupies the left (— ^) or rig ht (+|) V site of a rung and a denotes the z-component 
of the spin. Using these operators we define the conditional creation operators 



«L = (1 - ^ist) (1 - nio^i) (1 - rii^^) a\^„ 
similar to the operators used in the t-J-model PSI. We then obtain for Hpp: 



Hpp = Hs + J2 2^^^ + E ^il't^T^ 



where the K, Hs and tp contain spin operator products. It turns out to be useful to work 
with spin and pseudospin operators 

(T1,(J2,CK 



2 

CT,ai,«2 



In equation Si denotes the spin of the electron on the i-th rung and the z component of 
Tj corresponds to a, so the pseudospin Tj describes the hopping of the electron between the 
two V sites of a rung. Furthermore 
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4SiS, - 1 



ASiSj - 1 



pp\ 



2{Vn + Al^-El 
U + AY,-Er Vn + Aj.-Er 



2tl [aS.,S, - l) 2tl [aS^S, - l) 



U + A^j^-EPP U + Al-EPP 

D 



Vr + AI-E^ Vn + Al-Er 
^ . . 2tl (^4S,S, + l) 2tl (iS.S, + l) 
^^'^ = \/k + Ar - Er ^ Vn + Al- Er 
_ . 2tl Us J, + l) 2t|, USiSj + l) 

— £ !^ L (A) 

"^^y Vn + Al-Er Vn + Al-Er' 

Here A^ and A^ are additional energy differences due to the local change of occupation 
numbers with hoppings. Klz and Kil describe the most important interactions between 
electrons on neighbouring V sites. They drive the system to local "zig-zag" type ordering or 
local "in-line" type ordering of the V 3d electrons. Assuming that we have locally complete 
"zig-zag" short-range "ordering" within a ladder and no correlations between ladders the 
A^, A^ are given by 



(5) 



Assuming that we have local "in-line" "ordering" , they are equal to 

Al = ViL - Vl, Al = 2ViL - 2Vl, 

AY, = VjL - 2Vl, Al^ = 0, 

A^ = -2Vj,, Ag = VjL, 

A7^ = -ViL, A^^ = 0. 





= Vl, 


A]^ 


= 0, 




= m, 


AId 


= 0, 


K 


= -Vl, 




= 0, 


Ay. 


= -ViL, 


AL 


= -Vl. 







We now make a mean-field approximation for the spin operator products contained in 
K and tR. With this we obtain an effective pseudospin model for the charge degrees of 
freedom. K = {K) and = {tfj} become effective pseudospin coupling constants. This 
approximation should be acceptable since the corrections of the pseudospin part due to 
spin-spin-interactions are moderate. 
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With this approximation Hpp becomes the Hamiltonian for a strongly anisotropic pseu- 
dospin Heisenberg model in an "external field" 2t/j. We calculate the effective coefficients 
t/j, Kl, Kjl assuming that (S'jS'j+i) = — | along the ladders and (SiSj) = | between lad- 
ders. This assumption is based on the LDA+U results for the sign of the intra-ladder and 
inter-ladder exchange constants We use A^, A'^ from for the case of local "zig-zag" 
"ordering" and set E^^ to be the ground-state energy for PHP. Eq^ is then given by the 
expression for the ground-state energy of the Ising chain in a transverse field represented 
by tR. 



Eo = 2 i2tp) J '-^E 2J-—^ (7) 




with Ao = VL/2tp and E the complete elliptic integral of second kind. From the energies 
calculated in p6| by LDA-I-U for different configurations one obtains 



2Vl - VjL = 0.027 eV. (8) 
To find Vfi we assume that it can be obtained from Vl by scaling by with the different 



tR = 


-0.190 eV 


KiL = 


-0.677 eV 


Kl. = 


0.679 eV 


Kl. = 


-0.027 eV 


KLy = 


0.011 eV. 



distances d of the V sites. Using the values for the parameters from Ref. |T^, i.e., tji = 
-0.172 eV, ti = -0.049 eV, to = -0.062 eV, tjL = 0.110 eV, Vl = 2tR and f/ = 4 eV, we 
find 



(9) 



and Kiy are smaller than the other parameters by more than one order of magnitude 
and will therefore be neglected. Thus our model consists of Ising interactions Klz, Kjl and 
a "transverse field" 2tR 

Hi= Yl KlzTIT^+ KiLT:T; + Y.^tRT:. (10) 

<«.i>i {i,j)iL i 

The resulting geometry is triangular, as can be seen in Fig. |^a. Note that the first term is 
antiferromagnetic and the second term is ferromagnetic, resulting in geometrical frustration 
for the first two terms of equation (|). 

We first consider the case tp = 0. For this case the the relation between the absolute 
sizes of Klz and Kjl determines the behaviour of the system ||3^. If \Klz\ < \Kil\ the 



system undergoes a phase transition at some finite temperature into a low-temperature 
"ferromagnetic" state of the pseudospins with 2D-long range order. This state corresponds 
to an "in-line" -ordering of the electrons, i.e., chains of V^"*" and V^^ alternate along a- 
direction. 

If \Klz\ > \Kil\ the system remains disordered at any finite temperature. At T = it 
enters an antiferromagnetic state of the pseudospins with ID-order along b-direction and no 
correlation between neighbouring ladders. This corresponds to a "zig-zag" -ordering of 
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and V^^ along the ladders. Thus the system is effectively one-dimensional, although there 
are correlations between next-nearest-neighbour ladders . 



In a'-NaV205, we have \Klz\ > Therefore there is no phase transition at finite 

temperatures within the model described by equation ([l0|), provided that f/j = as assumed 
here. This is consistent with our choice of A^, A*^. This qualitative result remains true for 
other choices for the hopping integrals resulting in different Klz, Kh, e.g. those obtained 
from LDA in PB|. For tn = the system is therefore close to a quantum critical point with 



a transition from ID ordering in b-direction to 2D long range order \p2^ . 

Next we consider the case tji ^ 0. We note that the total Hamiltonian Hi does not imply 
a phase transition into an ordered state at any finite temperature. This is due to the fact 
that the first two terms of Hi do not lead to a phase transition as argued above and that the 
nonzero transverse field suppresses ordering even more. This result remains qualitatively 
true if we include the Kl^ and K^y terms from (P) which also suppress ordering. It also 
remains true if we include a small interlayer Ising-like interaction. This demonstrates that 
the phase transitions observed in a'-NaV205 cannot be of a purely Coulombic origin: to 
understand what happens in this material at low temperatures we have to include at least 
the spin degrees of freedom beyond the mean-field approximation used above. Furthermore, 
we know from Raman spectroscopy that there is a strong coupling between the charge 



distribution on the V sites and the distortion of the system |^3[. A change in the position of 
a V sites implies a change in the positions of the O site on the leg of the neighbouring ladder. 
The displacement of the O ions can in turn affect the effective superexchange interaction as 
observed from experimental data in . 



Therefore we generalize the Hamiltonian from equation (j^), such that it contains lattice 
degrees of freedom. Such a Hamiltonian has the following form: 

Hissp = ^ ^Lz'^ij'^i+ij + 2tR7^^ + ^ JijSijSi+ij 

+9ls J2 '^i/oUij + ^gf^b\,o^gjbvO,gj 
ij <l,j 

+ E <%a,jNa,g, + Yl ^>LlM^'^r (11) 



1,3 1,3 



Here i numbers the rungs on a ladder, i.e., it denotes the (pseudo)spin on a chain, j denotes 
the ladder/chain in the lattice and eo is a unit vector. We have introduced phonon operators 
b]jj, bgj for the displacement of neighboring V and O ions (VO) and for the displacement of 
Na or rung O ions {Or), u and Jij are given by 

J = f 1 + ^V-O ^VO ^ ^Na ^^Na ^On ^On 

^■J {Tij,Ti^ij) 

Xq K^N) ^ ^ 

Q.iXq) =h\^+h,,. (12) 
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Hissp essentially describes a spin-Peierls model with additional coupling to a one- 
dimensional Ising chain in a transverse field. It has to take into account the following 
observations. 

At T = K there is no correlation between pseudospins of neighbouring ladders for the 
pure Ising interaction case. Regarding the Coulomb interactions we may therefore treat the 
ladders as independent, setting Kjl = 0. The pseudospin part of the Hamiltonian (|^) for 
each ladder is then described by a one-dimensional Ising model in a transverse field yielding 
the first two terms of Hjssp- The property that the Hamiltonian (p!OD describes a model 
close to a quantum critical point has to be included in Hjssp- Therefore we have to use 
an effective value for Kl^ which is close to the value needed for quantum critical behaviour 
of the pseudospin Ising chains. This effective value Kl^ is not necessarily identical to the 
"true" value of this Coulombic coupling in the material. 

The spin part of the system behaves like a one-dimensional system at high tempera- 
tures i.e., we can use a one-dimensional Heisenberg model for each ladder to describe 
it. These terms are contained in the third term of Hjssp- The magnon dispersion at low- 
temperatures in a-direction is much smaller than in b-direction [Q,^ , even though interladder 
spin-spin coupling has a significant value p6|. Close to the disordered phase, however, in- 
terladder coupling is reduced due to frustration: as far as interladder spin-spin-coupling is 
concerned the system consists of triangles with one large antiferromagnetic and two smaller 
ferromagnetic interactions. 

On each ladder spin and pseudospin degrees of freedom are coupled. This coupling is 
incorporated into the effective pseudospin coupling Klz and the effective coupling constant 
Jij for the superexchange along the ladder. 

A shift of V and neighbouring O ions is assumed to cause an effective staggered field 
in z-direction for the pseudospins on one leg of a ladder described by the fourth term of 
Hissp- Such a shift also causes a change in the superexchange for the neighbouring leg of 
the neighbouring ladder (see Fig. ^ which has to be incorporated into the effective coupling 
constant Jij. 



As has been argued in [|l9l, a shift of the Na ions in c direction alternating along the 
ladder direction will cause an alternation in the superexchange along the ladder. The same 
effect is obtained by an alternating shift of O ions on a rung in b-direction. Such terms are 
included in Jij and their significance will be examined in section [111 C| . 

The last three terms of Hjssp describe the energies of the important lattice modes. 

The spin-Peierls coupling contained in the expression for Jij in (O) is slightly different 
from the normal form since it is the shift of the intermediate O ions or Na ions instead 
of the V ions, which causes the superexchange dimerisation. For parameters we will use 
tji = —0.190 eV and Jo = 0.045 eV for the case without charge ordering and no distortion. 

The effective superexchange along a ladder Jij defined in (|l^) consists of three parts: Jq 
denotes the original superexchange for the ladder without charge ordering or distortions. 

The second part, Jq f (T^^-, T^!,,^) describes the infiuence of the charge ordering on the 
superexchange. From numerical calculations we know that homogenous charge ordering 
on a ladder j causes a drop of the effective superexchange which is quadratic in (Tj), as 
can be seen in Fig. |^. The parabolic form of the curve can be understood from noting 
that the spectral weight of pair states with two particles on the same rung or the same 
leg decreases with 1 — (1 + 5co)(l ~ ^co) = ^co- Analytic calculations show a similar 
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behaviour [|3^. Therefore Jq f (Tfj,Tl^y) is approximately proportional to (Tj)^ and the 



square of the distortion accompanying the charge ordering. If we use the experimental value 
for the phonon frequencies this effect is already contained in the high-temperature values 
of the phonon dispersion Ugj and the term Jq f (T^^-,T^^y) is omitted, essentially making a 
mean-field approximation with regard to the spin operator products for this expression. 
Finally, there is the essential contribution 

(i + «Ti,vr4. + sfA.v-.,, + sf4/f4,) 4' (1 + / (Ti' 

to the superexchange. It describes the coupling between distortion, superexchange dimeri- 
sation and charge ordering. Here we use a mean-field approximation for the pseudospin 
operator products. The term causes a competition between the spin-Peierls order parameter 
and charge order parameter on a ladder. Klz will be fitted to properly describe the system. 
As argued above, we must use Klz instead of Kl^ in equation (111]). 



III. PHASE TRANSITIONS 

In this section we analyse the Hamiltonian given by equation ([11]). We have two types of 
instabilities to consider based on the third and fourth term. The first type results from the 
fourth term of equation ([TI|) which describes a coupling between local distortions and charge 
ordering. By distortion the system reduces the interaction energy. Due to elastic coupling 
this is also accompanied by an alternating exchange coupling on neighbouring ladders. The 
second type results from the third term in equation (|ll|) which includes a coupling between 
lattice and spin degrees fo freedom leading to an instability of a spin-Peierls type. 

In the following subsection we separate the system into two sublattices and first treat 
each subsystem on its own. Using RPA we analyze in subsection [illB| the first type of 



instability accompanied by charge ordering and superexchange alternation. We fit the free 
parameters gjs and K^z as to obtain the proper value of Tci and the shift of Tci in a 
magnetic field. By analyzing the ground-state energy of different ordered states we find that 
at low temperatures half the ladders are charge-ordered and on the remaining half there is 
superexchange alternation. 



In subsection |III Q we examine the influence of charge ordering on the second type of 



instability. The distortions of the charge-ordered ladders change the coupling u^^i j^f^i jJo 
and therefore the spin-lattice interaction. It will be shown that this interaction increases 
and results in a second phase transition. As shown below, this explains the anomalous BCS 
ratio. 



A. Separation of the subsystems 

The Hamiltonian (^) describes a model consisting of two subsystems 1 and 2 shown in 
Fig. 03. Both subsystems contain Ising pseudospin chains alternating with Heisenberg spin 
chains. The Ising pseudospin chains describe the charge distribution on a ladder, and the 
Heisenberg spin chains describe the spin on a neighbouring ladder. The degrees of freedom 
on each rung {ij) are characterised by a pseudospin Tjj on rung {ij) which is contained in 
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one subsystem and a spin Sij contained in the other subsystem. The two subsystems 1 and 
2 describe the complete spin-charge dynamics of the Hamiltonian ([TT|). 

Neighbouring chains of the same subsystem are coupled by the lattice dynamics. A 
change of the superexchange along a spin chain is obtained, e.g., by a displacement of an O 
ion on a leg. Via the elastic coupling this causes a displacement of the neighbouring V site 
(see Fig. ^). Such a local distortion changes the chemical environment and therefore the on- 
site energy of the vanadium dxy-oihiidX. This change in energy corresponds to a longitudinal 
field for the pseudospin of that rung. We assume that this field is proportional to the 
displacement dy of the V site of that rung and write: 

h^T = '^gisdv- (13) 

The two subsystems are coupled due to the dependence of the coupling Jij on the charge 
distributions {Tfj) and {TIj^^^). J^j describes the coupling of spins Sij and Si^ij in one 
subsystem and this term is contained in one subsystem, while {T^-) and (Z^+ij) are operators 
used in the description of the other subsystem. The coupling between the charge-ordering 
evidence by nonzero {T^-) and (T^+ij) and the superexchange coupling J^j along the ladder 
is negative. With increasing charge ordering along a ladder the effective superexchange 
coupling Jij decreases as shown in Fig. ^ This may change, though, if we include the 
effect of lattice distortions accompanying the charge ordering. We then have to deal with a 
combination of effects on Jij gued above. 

We first consider equation (0) for each of the two subsystems separately and then discuss 
the combined ordered ground-state. We can make a Wigner- Jordan transformation of the 
spin operators of a subsystem to obtain spinless fermions describing the spin degrees of 
freedom. Let kp = 7r/26 be the Fermi momentum of these spinless fermions where b is the 
lattice constant in b direction. 

For a single subsystem, e.g., subsystem 1, alternation of the superexchange Jij along the 
Heisenberg chains can be obtained by shifting the Na ions alternatingly in c-direction or by 
shifting the O ions on the rungs in b-direction as shown in Fig. ^. Applying the result of 
Cross and Fisher we arrive at renormalized phonon frequencies for q = 2kF'- 

^lk„Oa = ^lk„Oa - 0-26 \2gon (2A:f)|'T-\ 

^lk,,Na = ^lk,,Na " 0-26 \gMa {2kF)\^ (14) 

The coupling of a shift of the Na ion to a superexchange alternation is smaller by a factor 
of 2, because a shift of the Na site in c-direction affects the superexchange only on one 
two-rung cluster whereas a shift of the rung O site in b-direction affects the superexchange 
on two neighbouring two-rung clusters. 

For a single subsystem, e.g, subsystem 1, there is also an instability towards superex- 
change alternation along the Heisenberg chains accompanied by charge ordering along the 
Ising chains. This corresponds to shifts induced by the lattice mode go with go = 2/ci? 
as shown in Fig. |^. Within this mode the V sites on the chains of the A ladder shift in 
c-direction. The direction of the shift alternates along the ladder, corresponding to a stag- 
gered longitudinal field for the pseudospins of subsystem 1. The leg O ions on the B ladders 
also shift in c-direction together with the neighbouring V sites of the A ladder. This causes 
an alternation of the superexchange along the B ladder, corresponding to an alternation 
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of Jij along the Heisenberg spin chains of subsystem 1. Applying the results of Cross and 
Fisher and using RPA for the effect of the pseudospin-lattice coupling one obtains 

ul = ujI - 0.26 \gvo (go) ? T'' - Aglx.o (T) (15) 

where the last term describes the influence of the Ising chains. Here Xq (T) is the suscep- 
tibility of the Ising chain due to a longitudinal field along z-direction. gyo describes the 
coupling of a shift of the V sites on the superexchange of the neighbouring ladder via a shift 
of the O site. 



B. The phase transition at Tci 

We know from experiment that the phase transition at Tci in a'-NaV205 is accompanied 



by charge ordering B. We therefore first consider equation (p!5|). We can take g for the go 



mode in Fig. ^ directly from the results in |19]: 



vo ^ = 7.884 ■ 109 K s-i/2 (i6) 



^10 V my dy 

with dy = 6.0768 pm being the square root of the average of the squared shifts of the V 
sites on the A ladder, 6^^"^ = 0.26 the exchange dimerisation for layer a, J = 522 K and 
mv the V ion mass. We also need to obtain Xq{T)- While this is very difficult for general 
values of g, it is facilitated considerably when g = go- In this case the effective field for the 
pseudospins due to the lattice distortion is staggered in b-direction. The Ising chain has an 
antiferromagnetic coupling Klz > along b-direction. Calculating XqoiT) means therefore 
to calculate the susceptibility of an antiferromagnetic Ising chain in a transverse field for the 
application of an infinitesimal staggered longitudinal field. This is equivalent to calculating 
the susceptibility of a ferromagnetic Ising model in a transverse field for the application 
of an infinitesimal constant longitudinal field. This susceptibility can be expressed via the 
pseudospin correlation function p(A,T): 

1 ^ 

X,o (T) = P jim - J2 Pk-.I (A> T) (17) 

where P|'._^.|(A, T) = {TfTj){X,T) and (3 = l/ksT is the inverse temperature. A = Kizl^tR 
is the ratio between the pseudospin interaction and the transverse field and determines the 
properties of the Ising chain. \ = oo corresponds to the Ising chain without a transverse 
field, A = to the paramagnetic limit. 

For an Ising chain in a transverse field p^(A, T) can be written in the form of a Toeplitz 



determinant For large n acccording to Szego's Theorem ||3^ we have 



p:(A,T) = ip(A,T)G"(A,T) (18) 



where expressions for P{X,T) and G{X,T) are found in [0. Some simple limits are given 



below. To obtain x<?o(^) assume that the asymptotic expression (11^) is correct for all n. 



With this assumption summation of the series ( p^Tf ) yields 



1^^,, i + G(A,r) 

X,o,. = -/3P(A,T)^-^^. (19) 
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The approximation is expected to be correct close to the quantum critical point, i.e., close 
to A = 1 at T = 0, since then the susceptibility and the correlation length diverge. However, 
we can obtain P(A,T) and G{X,T) for the case A = oo, i.e., the pure Ising chain, from 
Ref. [071 as 



Gx=Z = tanh {pinX) . 
When A = we are in the paramagnetic regime and find 

Gx=o = A(tanh[/?t^(l + A)])exp((27r/?tH)"'/'exp(-/5tfi(l-A))...) (21) 
= 0. 

With these expressions we can calculate Xqoi'^) the limits X = oo, and find 



and 



X,o,A=oo(T) = ^expQ/5i^^,) (22) 

^^"''=° = 2^ ^^""^ (^^ ^^^'^O • ^^^^ 

Both expressions (|2^) and (^) are equal to the exact solution in these limits which implies 
that this Xgo,A should be a reasonable approximation for all values of A. We can therefore in 
principle calculate Tci from equation (^) knowing gjg, Uq^ and A by simply setting Ug^^ = 0. 

We can find from measurements of the elastic constants . There a strong anomaly 
in the Cge mode of a'-NaV205 was observed at the phase transition temperature. This mode 
couples to a zig-zag like charge ordering and we therefore identify it with the go mode of the 
present calculations. Using the high temperature value for the sound velocity fee = 4200 m/s 
we estimate = irvQ^/b = 279 K, where b = 3.611 A is the lattice constant in b-direction 
for the undistorted lattice. 

When is known we obtain A and gjs by fitting them to Tci = 34 K and to the 
experimental value of the shift of Tci in a magnetic field. The latter is reduced from its 
standard spin-Peierls value due to the Ising part of the transition. This is seen as follows. 
Considering only the spin-Peierls part of the transition, applying a magnetic field lowers the 
transition temperature Tci{H = 0) without a magnetic field to a value T^f^. On the other 
hand, the Ising chain susceptibility Xqo,\iT) increases with decreasing temperature. This 
causes an increase of the critical temperature from T^j" to the observed phase transition 
temperature Tci{H) > T^f^ in a magnetic field. 



The expansion for Tci in a magnetic field H for a pure spin-Peierls system is [39 



Tci(H = 0) \kBTci{H^O) 
Due to the influence of the Ising part this expression is modified to 
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tS{"{0) V^B^ci (0), 

'^^^ci (0) = 0.8— — (25) 



in the present case with 



(Tci m = u^l - Aglx,o,x (Tci (H)) . (26) 



For small shifts this leads to 

TciiH)~Tci{0) 0.36 / i^bH 



with 



Tci (0) 1 + Tci (0) a (Tci (0)) yksTci (0) 



(27) 



a (T) = ^ ^"^^ , (28) 



From the experimental result for the shift of Tci = 34 K in a mag netic field in i,P,|0| 
obtain a (Tci) ~ 0.086 K~^. Fitting gjs and A to these two values we obtain 



we 



A = 0.99850 

gis = 1.688 K/(pm shift of V). (29) 

Note that this value for A is close to the value for the ratio Kjl/Klz ~ 0.997 from equation 
(^). This indicates that geometrical frustration plays an important role. To see how stable 
the solution is towards a change of the parameters, we evaluate the dependence of Tci on 
Klz, f^qo and gis'. 



dT, 



ci 



dTci 



dT 



ci 



dg 



1.06, 
-0.14, 

12 pm. (30) 



Is 



Tci is rather sensitive to small changes of Klz or A. This sensitivity is due to the fact, that 
(1 — A) -C 1. It is much less and the stability therefore much improved compared to the 
result presented in [^, where A ~ 0.99986 was suggested. 



Next we include the coupling between the two subsystems. For this we will use the mean- 
field values for (Tfj), {T^^y). They are zero at Tci and therefore the critical temperature Tci 
remains unchanged in this approximation. But a mean-field approximation is not reliable 
close to the critical point of the transition where fluctuations must be taken into account. 
We therefore consider the zero temperature free energy of different ordered states. For this 
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we need the dependence of the superexchange Jq along a ladder on the charge ordering 
order parameter (T?) on the same ladder. This is approximately 

^? = ^o^(l-4(T,0') (31) 

where we neglect a small contribution to Jfj^ from the diagonal hopping in the case of com- 
plete charge ordering on the ladder, i.e., for (TJ) = |. As argued at the end of section ||, Jf^ 



enters our equations in two places. First, it effectively changes ajg^ through Jq/ (T^pTl^y). 
Since we used Ug^ from experiment, this effect is already included. Second, J^j^ modifies the 
energy gain resulting from the spin degrees of freedom due to distortions and it is this effect 



we have to consider now. Using the expression for the magnetic energy gain from Ref. [12 
for the alternating Heisenberg chain we can write down the energy gain dF{T = 0) of the 
ground-state of the ordered system versus that of the ground-state of the disordered system. 
The dimerisation parameter of the Heisenberg chains S^j^'^"^^ with 6^j'^'^'^^ = g^^Qi^2/ J from 
equation (p!6[) is the order parameter for each of the two subsystems. Q = dy^Jniv/h is the 
canonical displacement. We find: 

dF{T = Q)=b [Ql + Ql) - EsP m (1 - 4 {T^f) 

-EsP m (1-4 {Tlf) - 2gu m)Qi + {T^^) (32) 



with 



1 2 



Esp{Qi,2) = 0.3134J 



4/3 



{TI2) = ^ (1 - (-4X7. (T = 0) gisQi,2)) ■ (33) 

Here we assumed that approaches | exponentially if a staggered parallel field gisQ is 
applied. In equation (^) Esp denotes the energy gain from the spin-Peierls distortions, while 
the terms proportional to denote the energy gain from charge ordering. The first term 
stands for the elastic energy of the distortion. Optimizing equation (^) in (5j G [0, 1] 

we find a minimum at S'^P = 0.023, 6q^ = or vice versa with dF = —0.80 K. Requiring 
= Sf^ = Sj, i.e., the equivalence of the two subsystems we find 6j = 0.009 with 
dF = —0.35 K. We have plotted in Fig. ^ dF for 6^ = df'' and 5^ = 0. In the previous 
calculations we used S'^f^ = 0.26 to obtain (7^*^, so the first result is more consistent with 
our calculation than the second one and it also has the lower energy. 

From the derivation of xis we can obtain the bare correlation length along the Ising 
chains as 

At Tqi we find that ^ 183 lattice constants. This large value of C,is might offer an 
explanation for the strong dependence of Tci on doping. By substituting Na with Ca or 
depleting the system of Na additional electrons/holes are introduced into the system. This 
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implies double/zero occupancy of rungs. They weaken correlations along the Ising chains. A 
rough estimate of the critical doping value at which the phase transition temperature Tci is 
suppressed is xc ~ 0.5%. We obtain it by assuming xc^is ~ 1. This value is in qualitative 
agreement with experiment. Substituting Na by Ca the transition disappears between 1% 
and 2.5% of Ca [^. For hole doping due to Na deficiency the transition disappears between 
2% and 3% deficiency [||. 



C. The phase transition at Tc2 

Next we consider the second phase transition observed. According to this transition 
is of the Ginzburg-Landau type. It opens a spin gap and creates a local distortion at the Na 
sites as observed from the changes of the quadrupolar electric field tensor. Charge ordering 
at the V sites starts before this transition sets in ||^. For these reasons and in order to explain 
the observation of eight inequivalent Na sites in ^^Na-NMR we suggested in that 



the opening of the spin gap in the system is due to an alternating shift of the Na ions along 
the charge-ordered A ladders as shown in Fig. A second possibility is an alternating 
shift of the rung O ions on the A ladders as shown in Fig. ^d. 

By using the hopping matrix elements obtained in [IT3 for the charge-ordered ladder in 



layer a we can calculate the effective superexchange J dependence on the charge ordering 
6co on the same ladder. As we did for the undistorted ladders (see Fig. H), the superex- 
change is defined by the singlet-triplet gap on a two-rung cluster. It is obtained via exact 
diagonalisation. We know from experiment that J ~ 440 K in the low-temperature phase. 
This value of J corresponds to an incomplete charge ordering Sqq = \{n+ — n_) ^ 0.32 on 
the A ladders in agreement with experiment P]. 

The hopping matrix elements t^, tn, change due to the distortion accompanying the 
charge ordering [1^. We parametrize these changes with t{5co = 0) — t{5co) oc 5iat oc 5co, 



and normalize the prefactors such that the hopping matrix elements from Ref. for the 
high temperature undistorted ladder are obtained at 5co = and those for the charge 
ordered ladder are obtained at 5co = ^co- 

Due to alternating shifts of the rung O ions or the Na ions, the effective hopping along 
the leg and the effective diagonal hopping t/j alternate also. The size of this is given 
by the parameter ^ = — t^^)/{t^^ + t^'^). If we know 5t^ and we can find 
the superexchange alternation parameter 5j by calculating the superexchange for the in- 



equivalent two-rung clusters |T^. This gives us an approximate value for the spin-lattice 
coupling. Using it, the parametrized t^, tjj and tji from above, and equation (0) for the 
critical temperature we find Tsp{Sco) for a spin-Peierls transition depending on the charge 
ordering: 



J., , - 0.8 ;7;7" . (35) 

J{6co) T^J[5co)^c 



We can estimate the value ujqx from the velocity V22 of the C22 mode in [^]. It describes a 
longitudinal mode along b-direction as required for the proposed shifts of the rung O ions 
with V22 = 6500 m/s which yields uqx = 432 K. For given values of and we obtain the 
dependence of the critical temperature Tsp on the charge ordering 6co as shown in Fig. |7]. 
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In the low-temperature limit the spin gap is approximately given by A{Sqq) ~ 1.8Tsp{Sqq). 
Therefore the BCS ratio is enhanced by a factor of Tsp{S^q)/Tci- 

We search for values of and for which the following conditions are fulfilled: (i) 
without charge ordering the critical temperature for such a transition is below 34 K, since 
the spin gap opens only below the charge ordering temperature; (ii) with increasing charge 
ordering the critical temperature is required to increase, because an enhanced BCS ratio 
has been observed; (iii) at 6co = ^co ^Pi^ S^P ^('^co) ~ l-8^c^('^co) should be 
approximately 100 K as found experimentally. 

The result is found in Fig. |^. The conditions (i) to (iii) are fulfilled for = —0.059, 
6tjy = 0.090 and with Tci = 34 K an effective BCS ratio of 6.05 is found. We therefore 
have a spin-Peierls phase transition which is driven by charge ordering on the same ladder. 
When the first transition takes place at about Tci = 34 K charge ordering sets in on the 
A ladders. This changes the two-particle correlation functions and the effective hoppings so 
that Tsp increases as compared with its value without charge ordering. At some temperature 
Tc2 < Tci we have Tsp{Sco(Tc2)) = Tc2 and a spin-Peierls transition driven by charge 
ordering takes place. This is caused by an alternating shift of the O ions on the rungs of 
the A ladders, probably together with a small alternating shift of the Na ions along the A 
ladders. This also causes an inequivalence of the Na sites along the A ladders, such that 
we have 8 inequivalent Na sites as observed in experiment Since this happens before 



6co attains its final value, the spin gap at lower temperatures is larger than what would be 
expected from the standard BCS ratio and Tci- 

In order to see whether the values for and are plausible, we have investigated the 
effects of alternating shifts of the rung O or the Na sites on the effective hopping matrix 
elements between V sites on the ladders. For this we use the Slater-Koster method as 
described in [|TU| . The results are given in Table |I[ Those for and have been obtained 
for otherwise undistorted ladders. and are found to have opposite signs, as assumed 
above. This can be understood by noticing that an increase of diagonal hopping due to 
increased hopping via the Na sites goes together with a decrease of the hopping along the 
legs of the ladder. The effects of charge ordering distortions and of the shifts of Na or rung 
O sites on the hopping matrix elements are small. Therefore we may neglect higher-order 
effects resulting from combinations of the distortions. 

The results in the first three rows of Table | are based on [|1^] where we included V- 
O, V-V and Na-0 hopping matrix elements in the initial Hamiltonian, projecting these on 
effective V-V hopping matrix elements. However, it has been argued in PB[ that direct 



hopping between the O sites should also contribute to the effective Hamiltonian. It remains 
unclear, which fraction of the total 0-0 hopping matrix elements obtained in |]2B| results 
from direct hopping and which fraction comes from indirect hopping via the neighbouring 
Na site. To investigate the effect of such a hopping on the shifts of the rung O ions we 
introduced 0-0 hopping matrix elements according to the Slater-Koster approximation for 
into the initial Hamiltonian and calculated effective V-V hopping matrix elements as before. 
One finds tp = -0.17 eV, = 0.17 eV, tn = -0.04 eV, tji = -0.16 eV. Introducing an 
alternating shift of the rung O ions in b-direction by 1 pm leads to ~ 0.003 and 0.1 
as given in the fourth row of Table |I[ The effect especially on the diagonal hopping is much 
stronger than it is without inclusion of the direct 0-0 hopping matrix elements. However, 



these values for and tj^ are far from those found by LDA [BH] and ab-initio methods [HS 
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One can therefore conclude that at least within the Slater-Koster method the contribution 
of the direct 0-0 hopping should be small as expected for next-nearest neighbour hopping, 
although it may help to explain the small differences between the LDA results of |^ and 
our results in The results from Table | can serve only as an estimate of the values of 
the spin-lattice coupling. 



IV. COMPARISON WITH X-RAY DIFFRACTION AND NEUTRON 
SCATTERING EXPERIMENTS 

In the following we compare our results with experimental observations. We do this 



with respect to two types of experiment: x-ray structure determination |[T7HT9| and inelastic 
neutron scattering results for the magnon dispersion compute the latter for the 

calculated structure by using a local-dimer approximation. 

Comparing our results with x-ray structure determination we find good agreement. In 
both cases we have a superstructure consisting of two inequivalent ladder types alternating 
along a-direction: charge-ordered A ladders and strongly dimerised B ladders. A phase with 
two A ladders separated by a B ladder is also correctly found. We can explain these findings 



with the analysis given in section |111 B|: the displacement of the V sites on the A ladders and 



of the O sites on the B ladders cause a combined charge-ordering-spin-Peierls transition. 

X-ray structure analysis finds Fmm2 symmetry in the low-temperature phase. In our 
model this symmetry is realized for Tci > T > Tc2- However, the phase transition at Tc2 



obtained for our model in section |III Q is accompanied by a breaking of this symmetry. It 
results from the shift of the rung O ions on the A ladders, as shown in Fig. A low- 
temperature spin gap of A «i 100 K on these ladders requires a superexchange dimerisation 
of about 0.05 when A = 2J5^J^ from Ref. [|1^] is used with the experimental value for J = 440 
K. The coupling constant at Sqq = 0.32 is approximately 0.043 per pm shift of the rung O 
ion. This corresponds to an actual shift of 1.1 pm of each rung O ion on ladder A. 

The intensity of the x-ray scattering for small scattering angles is proportional to 
where Z is the ionic charge. Below Tc2 only one O ion per two formula units shifts its 
position to break the Fmm2 symmetry. In addition as estimated above these displacements 
are smaller than the displacements caused by the phase transition at Tci and below, where a 
displacement of 7.46 pm and 4.26 pm for the V sites on the A ladder is found experimentally. 
It is therefore well possible that the scattering peaks resulting from the lower symmetry 
cannot be observed within experimental resolution. Insofar x-ray structure determination 
results do not contradict our theoretical results for the second transition. 

Magnon dispersion, however, should be greatly affected by a dimerisation. The structure 
which is obtained after the two transitions have taken place is schematically shown in Fig. p. 
The calculated magnon dispersion for this structure should therefore be a good test of 
the theory. For other suggested structures it was found that the corresponding magnon 
dispersion disagrees with experimental results: in the case of a lattice with "zig-zag" charge 



order on all ladders, the dispersion along a-direction does not agree with experiments ffi 



For the spin-cluster model proposed in [T^ the magnon dispersions, both in a- and b- 



direction, do not agree with experiments either P5|-^]. 

We therefore calculate the magnon dispersion at a temperature well below Tc2- In 
order to do this we use the spin-dimer representation following Refs. The initial 
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magnetic Hamiltonian is then 



+ 5^JB(l + (-iy(5B)4-4+l 



{ij,mn) 



for a geometry and J)-^"^" as shown in Fig. It describes the pure spin part of equation (|TT]) 
with the effects of lattice distortions and pseudospins included within the effective coupling 
constants. The large value of 6b ~ 0.26 at T = 15 K found in corresponds to a large 



gap in the excitation spectrum of the B ladders. Using = 2Jb6^J'^ for the spin gap |T2 



and Jb = 52 meV from Ref. [TU[ we find fti 38 meV. The B ladders should therefore 
not contribute directly to the low-lying parts of the magnon dispersion. Since we are only 
interested in these, we assume that we have an indirect exchange coupling between A ladders 
through the virtual singlet-triplet excitations of the B ladders. This results in an effective 
Hamiltonian: 



H — y ^ Ja ^1 + (~1)'' 6a) Sij^Sij^i^ + y ^ JaSij^Si^ij^ 

ij-y ij-y 

+ y ] Jr)Sij-y ^Sj+ij+i^ + Si^ij-i^) 

~l~ ^ ^ JcSjj-ySjj^ (37) 

with corresponding geometry and couplings as shown in Fig. This exchange Hamiltonian 
exhibits manifestly a doubling of the period along the a axis. Here we also introduced a 
superexchange Jc in c-direction between layers; it is assumed to exist mainly between the 
V21 sites which lie directly above each other [T^. We therefore have 7 = ±1: Along c- 
direction the V sites of the A ladder are ordered as . . . -V21-V21-V22-V22- ... as shown 
in Fig. |lO| |ir7|-[l9(|. Next we introduce dimer variables. We denote each dimer within the A 



ladder (see Fig. ^ by coordinates ij'y, which are not identical with the former coordinates 
of the spins. As denoted in Fig. ^ we then have spins Sij^i^2- The dimer variables are then: 

Kij-y = Sjj^l + Sij^2 ^gg-^ 
Lij'y Sij-yl 

With those, equation (^) takes the form 
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H — — J (^1 -\- 6a) ^ ^ Kij^Kij^ Lij^Lijy 

ij-y 

-\- — Ja (1 ~ 6a) ^ ] Kij^Kijj^i^ — Lij^Lij^i^ 



4 



«J7 



1 >^ - - 



+-^i+lj+l7 ~ -^i+lj+l7^ + ^-^ij7 + -^ji7 



2 



As in Ref. we will only use the products Lij^LmnY, since the other terms do not con- 
tribute to the dispersion of the spin excitations. We then transform the L to 



4-7 = 75 ( + 
47 = 71 (^-^ 



(40) 



Furthermore we assume that the ladder designated by (i + 1,7) is equivalent to the ladder 
(i,7). This corresponds to a geometry of the A ladders shown in Fig. We therefore set 
Lij^ = Li^ij^, and separate the Hamiltonian into two parts: 

Hm+ = -\ {J A (1 + Sa) - Jc) MtjMtj 

-\ja (1 - 5a) Mt,M±^, - haMt^Mli, 



1 



Hm- = {Ja (1 + Sa) + Jc) M-^M-^ 

-\ja (1 - 6a) M-,M-,+i + \jaMr,Mr^,j 
1 



-Jo [M-^Mr^,^ + M^Mr^,^^,) . (41) 

The first term of each Hamiltonian describes an effective dimer with interaction strength 
Ja {1 + Sa) Jc- Its dynamical susceptibility is 

{Ja (1 + Oa) =F Jc) - 

The dynamical RPA susceptibility of coupled dimers without anisotropy of the superex- 
change is given by 

X (q, uj) = [l- J{q) u^x (^^)]~^ Uxx (t^) (43) 
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where J is the exchange term between dimers as given by the last four terms of and 
Hm-- After Fourier transformation these, for Hm+ and J~ for Hm- are given by 



{q) = ^Ja (1 - 5a) cos 2ky 



(44) 



Here kx = 2'n'qx/a, ky = 2Tiqy/h with a and h being the lattice constants of the high tem- 
perature unit cell. We note a doubling of the unit cell in b-direction as required before 
by equation (|37|) . The spin- wave excitations are obtained from the poles of x(g', cj). With 
equations (^2]) and (H3) we obtain 



Ul = {J A (1 + Sa) T Jcf - {J A (1 + 5a) T Jc) 
X {J A (1 — 5 a) cos 2ky ± {Ja COS (^2,. — ky) 

— 4: Jd COS kx COS ky)) . (45) 

We set J^^ = Ja — 4:Jd and J4 = 440 K, as used in our previous analysis and obtained from 
experiment 0. Then we fit the unknown parameters J^^ , Jc and 6a to the experimental 
values for three of the four gaps obtained from u;+ and u;_. These gap values have been 
observed for T < 4.2 K at g^; = 0, = | and qx = 0, qy = 1 0. Fitting three of the four 
gaps of sizes 10.9 meV, 10.1 meV, and 9.1 meV we obtain 



Jf = 0.212 meV, 

= 0.431 meV, (46) 
5a = 0.0310. 

This determines the fourth gap to be 8.14 meV which is in excellent agreement with the 
experimental value 8.2 meV from Ref. 0. In the model presented here these spin gaps 
result from the assumed dimerisation of the charge ordered A ladders whereas previously 
the origin of the spin gaps was suggested to lie in an anisotropy of the superexchange . 



The dispersion of cu^ along a-direction for g^^ = |, 1 are shown in Fig. |TT|. These curves also 
agree with experiments. 

The low value of J^^, which causes the dispersion along a-direction is due to an effective 
coupling of next-nearest-neighbour ladders. Furthermore, Ja is partly compensated due to 
Jd couplings. 

The value for Jc, which causes the gap between the two modes, corresponds to a hopping 
tc ^ 0.021 eV between neighbouring V sites of neighbouring layers when J^ = ^ and U = 4 



eV are used. This is in reasonable agreement with the value tc = 0.015 eV found in ||T9 
from the low-temperature structural data by a Slater-Koster approximation. 

We can also obtain the approximate dispersion in 6-direction for g^. = as shown in 



Fig. |T^. The maximum is at ~ 55 meV and agrees nicely with the value Umax = 59.5 meV 
estimated in from experimental data. 

We therefore conclude that the results of our analysis in section |ITT| agree well with both 
the experimental results of x-ray structure determination and the magnon dispersion found 
by inelastic neutron scattering results for the magnon dispersion. 
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V. SUMMARY AND CONCLUSIONS 



In this article we provided a theoretical description of the phase transitions in a'-NaV205. 
In section |I| we started from a single band extended Hubbard model describing the hopping 
of electrons between V sites. We projected this model onto an effective spin-pseudospin 
Hamiltonian similar to that in Refs. [p^ -p5|. Using parameters obtained from a previous 
Slater-Koster analysis and findings of LDA+U we argued that within the given model a 
phase transition cannot occur from Coulomb interactions only: the effective lattice for such 
a transition is triangular and its geometric frustration suppresses charge order, a distinctive 
feature which is present in the observed transition. We therefore constructed a minimal 
microscopic model based on a spin-pseudospin Hamiltonian which incorporates the coupling 
between charges, spins and the lattice. 

In section |T| we analyzed this model with regard to phase transitions. We divided the 
model into two subsystems. Each of it describes the charge degrees of freedom on the rungs 
of half the ladders and the spin degrees of freedom on the other ladders. For each subsystem 
we find an instability towards a transition at Tci which causes "zig-zag" charge ordering and 
superexchange dimerisation. With this combined spin-Peierls-Ising transition we can explain 
the anomalous shift of Ta in a magnetic field. It is reduced from its normal spin-Peierls value 
due to the influence of the charge ordering. By calculating the free energy at T = we found 
that the system should enter a phase where half of the ladders are "zig-zag" charge-ordered 
(A ladders), while the other half (B ladders) shows a strong superexchange alternation. The 
reason for this asymmetry lies in a competition between the order parameters for charge 
ordering and superexchange dimerisation due to a shift of the leg O sites on the same ladder. 
This explains observations of x-ray structure determination. Due to the strong effect of the 
one-dimensional Ising chains we can also explain qualitatively the strong dependence of the 
transition on depletion of Na or substitution of Na with Ca. 

By analyzing the coupling between superexchange dimerisation and lattice distortion in 
the charge-ordered ladders we found that another transition is induced. Charge ordering 
increases the coupling between an alternating shift of the O sites on a rung and the superex- 
change alternation beyond the threshold value of a second transition at Tc2- The system 
then enters a phase where the charge-ordered A ladders dimerise and a spin gap opens. Due 
to the alternation the Na sites along the A ladders become inequivalent, such that we have 
8 inequivalent Na sites as observed in ^^Na-NMR |]5|,|42[|, not only 6 as would be the case in 



the Fmm2 symmetry indicated by x-ray structure determination. 

This also helps to explain the discrepancy between measurements of the critical exponent 
of the lattice distortion [P-[TT1| and of the critical exponent of the spin gap opening the 
former is smaller than the latter whereas one would expect the opposite for a single transition 
due to A oc Sfl^. From the observations of a logarithmic peak in the specific heat at Tci and 
the observation of fluctuations by x-ray diffuse scattering [|TT| we conclude that the transition 
at Tci is mainly of 2D Ising character whereas the second transition at Tc2, which opens 
the spin gap, can be described within a mean-field theory. Since the charge ordering is 
not complete, when the second phase transition is triggered, the coupling constant further 
increases for decreasing temperatures. This leads to an additional increase of the spin gap 
and also of the BCS ratio. The conventional BCS ratio does not account for the temperature 
dependence of the coupling constant g^^^. 
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Within experimental resolution we find agreement with the x-ray structure determina- 
tion. The second transition breaking the Fmm2 symmetry comprises only a small shift 
of every tenth O site. We also find excellent agreement with experimental results for the 
magnon dispersion at a temperature well below Tc2- 

In conclusion we presented a microscopic model for Q;'-NaV205 which yields two phase 
transitions close to each other. The first is a "spin-Peierls-Ising" transition caiising charge 
order and superexchange dimerisation. The second is a pure spin-Peierls transition trig- 
gered by an increase of the coupling constant due to charge ordering. With this model we 
can explain qualitatively and quantitatively a number of experimental observations, i.e., the 
existence of two transitions, the general structure of the low-temperature phase, the anoma- 
lous shift of Tci in a magnetic field, the anomalous BCS ratio, the strong dependence of the 
transition on doping, and the observed low-energy magnon dispersion. 
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TABLES 





r tl3 + ~tD- 

tD + +tD- 




K J+-J- 


Na shift in 








c-direction 


0.0105 


-0.0135 


0.0035 


V shift in 








b-direction 


0.0236 


0.0136 


0.0347 


Rung shift 








in b-direction 


0.0072 


-0.0093 


-0.0016 


Rung shift 








in b-direction with 








too included 


0.102 


0.003 


0.096 



TABLE L Change of hopping matrix elements for an alternating Na shift in c-direction, an 
alternating V shift in b-direction and an alternating shift of rung O ions in b-direction respectively. 
Shift is by 1 pm per site on an otherwise undistorted ladder. The accompanying change of the 
superexchange has been calculated as the singlet-triplet gap on a cluster of two neighbouring rungs 
of a ladder p|. 
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FIGURES 




FIG. 1. Hopping matrix elements and Coulomb interactions in the V layers 
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FIG. 2. a) Lattice geometry if one introduces Sij and Tij operators and represents each rung by 
a point, b) The two subsystems of the model described by (^Tj) consisting of Ising (IS) pseudospin 
and Heisenberg (HB) spin chains. 
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FIG. 3. a) Alternating Na shift in c-direction schematically shown by arrows causes superex- 
change alternation, b) Alternating shift of the rung O ions in b-direction causing superexchange 
alternation and inequivalence of the Na sites along the ladder. 
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FIG. 4. The qq mode is shown, for which the V ion on one leg and O ion on neighbouring leg 
move together in c-direction as shown by arrows. This causes charge ordering on A ladders due to 
a staggered change of the V on-site energies and superexchange alternation on B ladders due to O 



ion motion. In experiment [19| shifts in a-direction are also observed which are not shown here. 
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FIG. 5. Dependence of superexchange J on charge ordering. Results from cluster calculation 
for two rungs of the same ladder with parameters tn = 0.172 eV, ti = 0.049 eV, to = 0.062 eV, 
Vl = 0.344 eV, Vr = 0.398 eV, ?7 = 4 eV and superexchange defined as the singlet-triplet gap |19|. 
"Zig-zag" charge ordering or "in-line" charge ordering was induced by changing the on-site energies 
with ei = £3 = —62 = —64 for "in-line" ordering and ei = —€2 = —£3 = £4 for "zig-zag" ordering. 
Inset shows geometry of cluster. The difference in the behaviour of the superexchange between the 
two ordering patterns comes from the fact that / t_D and Vl i^Vo = 0. 
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FIG. 6. Ground-state energy gain dF of the ordered system versus the disordered system 
depending on the order parameters. 
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25^Q=n^-n 



FIG. 7. T^p-dependence on charge ordering for a spin-Peierls transition due to alternating 
shifts of the rung O ions on the A ladders in b direction. Dimerisation of hopping matrix elements 
for a shift of 1 pm given by = 0.090 and (5j^ = —0.059. x mark conditions (i) and (iii) (see text). 
The increase of Ts p for small 5co is due to increase of the effective 6j resulting from increase of the 
hopping matrix elements due to the distortion accompanying the charge ordering. The decrease 
for large 6co is caused by the decrease of the average superexchange J as shown in Fig. 




A B A B 

FIG. 8. Geometry of initial magnetic Hamiltonian (^). J± = J(l ± 5) for the respective ladder. 
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A B A B 

FIG. 9. Geometry of effective magnetic Hamiltonian (|37|). 1 and 2 denote the two spins of 
a dimer, and Ju denote the effective superexchange between spins in the A ladder via the B 
ladder for nearest interladder neighbours and next-nearest interladder neighbours. V21 sites with 
increased electron density denoted by black circles, V22 sites with decreased electron density by 
white circles. 



30 




B A B A 



FIG. 10. Unit cell of low-temperature ordered state. Along the A ladders one has pairs of V21 
sites (black circles) alternating with pairs of V22 sites (white circles) in c-dircction. J^. therefore 
connects two layers 7, 7 only. For such two layers, dimers on the A ladders are assumed to lie 
directly above each other in c-direction. 
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FIG. 11. Magnon dispersion along a-direction for qi, = ^ = Q^^ (solid lines) and qb = I = Q 
(dashed lines). 
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